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The Fermi Large Area Telescope (LAT) has provided the measurement of the high energy cosmic 
ray electrons plus positrons (CRE) spectrum with unprecedented accuracy form 20 GeV to 1 TeV. 
Recently this range has been extended down to ~ 7 GeV. The spectrum shows no prominent 
features and it is significantly harder than that inferred from several previous experiments. While 
the reported Fermi-LAT data alone may be interpreted in terms of a single (electron dominated) 
Galactic component, when combined with other complementary experimental results, specifically the 
CRE spectrum measured by H.E.S.S., and especially the positron fraction measured by PAMELA, 
an additional electron and positron component seems to be required. We show that the acceleration 
of electron-positron pairs in Galactic pulsars may offer a natural (though not unique) interpretation 
of all those results. 



I. INTRODUCTION 

Prior to 2008, the high energy electron spectrum 
was measured by balloon-born experiments [l[ and 
by a single space mission AMS-01 ^2^. Those data 
are compatible with a featureless power law spectrum 
within their errors. This is in agreement with theoret- 
ical predictions (for a recent review see (sj) assuming: 
i) that the source term of CR electrons is treated as a 
time-independent and smooth function of the position 
in the Galaxy, and the energy dependence is assumed 
to be a power law; ii) that the propagation is described 
by a diffusion-loss equation whose effect is to steepen 
the spectral slope respect to the injection. Possible 
deviations from a simple power law spectrum may, 
however, be expected above several hundred GeV as 
a consequence of synchrotron radiation and Inverse 
Compton (IC) energy losses which, at those high en- 
ergies, limit the electron propagation length to a dis- 
tance comparable to the mean distance between astro- 
physical sources 0, @ or because the possible presence 
of exotic sources. 

Almost a year ago, the ATIC balloon experiment 
found a prominent spectral feature at around 600 
GeV in the total electron spectrum. Furthermore, the 
H.E.S.S. 0, [1] atmospheric Cherenkov telescope re- 
ported a significant steepening of the electron plus 
diffuse photon spectrum above 600 GeV. Another in- 
dependent indication of the presence of a possible de- 
viation from the standard picture came from the re- 
cent measurements of the positron to electron frac- 
tion, e+ /(e^+ e+), between 1.5 and 100 GeV by the 
PAMELA satellite experiment fgj. PAMELA found 
that the positron fraction changes slope at around 10 
GeV and begins to increase steadily up to 100 GeV. 
This behavior is very different from that predicted for 
secondary positrons produced in the collision of CR 
nuclides with the interstellar medium (ISM). 

The experimental information available on the CRE 



spectrum has been drastically expanded as the Fermi 
Collaboration has reported a high precision measure- 
ment of the electron spectrum from 20 GeV to 1 TeV 
performed with its Large Area Telescope (LAT) [l3] 
based on the first 6 months of data taking. Prelim- 
inary updated Fermi-LAT measurements, based on 
the first 12 month data, basically confirm the high 
energy spectrum measured by Fermi above 20 GeV 
and extend it down to almost 7 GeV. In that low en- 
ergy range the Fermi-LAT data are compatible with 
those of AMS-01. The electron -I- positron spectrum 
measured by Fermi-LAT is compatible with a simple 
power law of slope 7 = —3.08 and presents some hints 
of a hardening at around 70 GeV and steepening above 
~ 500 GeV. Although the significance of those fea- 
tures is low within current systematics, they suggest 
the presence of more components in the electron high 
energy spectrum. It is also worth noticing here that 
the hard electron spectrum observed by this exper- 
iment exacerbates the discrepancy between the pre- 
dictions of standard CR theoretical models and the 
positron faction excess measured, most conclusively, 
by PAMELA This makes the exploration of some 
non-standard interpretations more compelling. 



II. CONVENTIONAL INTERPRETATION 

In this section we discuss a simple possible inter- 
pretation of Fermi-LAT CRE data in terms of a con- 
ventional [2^ Galactic CRE scenario assuming that 
electrons sources are continuously distributed in the 
Galactic disk and that positrons are only produced 
by the collision of primary CR nuclides with the in- 
terstellar gas. To this purpose we use the GALPROP 
numerical CR propagation code [ll| . We start consid- 
ering a conventional model assuming a single power- 
law dependence of the diffusion coefficient on rigidity 
with index 5 — 0.33 (Kolmogorov like) and electron 
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(e^ only) injection spectral index 70 = —1.6/ — 2.42 
below/above 4 GeV. In [m] this model was shown to 
provide a good description of published Fermi-LAT 
data above 20 GeV (few other combinations of S and 
7o were also shown to fit the data). 

That model, however, faces a series of problems 
when compared with other experimental data and, 
most seriously, with the preliminary low energy Fermi- 
LAT, namely: i) below 20 GeV that model predicts an 
excess of electrons + positrons respect to Fermi-LAT 
and AMS-Olf^ data when a reasonable solar modula- 
tion is considered (see red line in Fig. [T]); ii) the model 
overshoot H.E.S.S. data above 1 TeV; iii) the positron 
fraction e^/(e+ -I- e~) they predict is not consistent 
with that measered by PAMELA Q (see the insert in 
Fig. [1]). 

While (ii) may be interpreted as a consequence of 
the stochastic nature of astrophysical sources (see Sec. 
2.2 in [3] and ref.s therein) or a high energy cut-off 
in the source spectrum, the other caveats are more 
serious and call for some revisions of the model. 

Disregarding high energy (E > 10 GeV) PAMELA 
data in this section, we have only to bother about 
the discrepancy of the conventional model with low 
energy data. It should be noted that below 10 GeV 
solar modulation is expected to affect significantly the 
spectra of electrons and positrons reaching the Earth. 
In the most simple scenario modulation is treated as a 
charge independent cooling effect which can be quanti- 
fied in terms of a parameter $ having the dimensions 
of rigidity (force field approximation). Such param- 
eter is time dependent and it is generally tuned to 
reproduce CR proton data. In periods of low solar ac- 
tivity (as those in which Fermi, PAMELA and AMS- 
01 have been operating) its value typically ranges in 
the interval 400 — 600 MV. None of those values, how- 
ever, allows the conventional model considered in [l3] 
to perfectly match Fermi-LAT data (see red lines in 
FiglD). Assuming a steeper injection spectrum below 
few GeV's, while allows to improve the fit of Fermi- 
LAT data, it overshoots considerably AMS-01. Such 
discrepancy may hard to be compensated even con- 
sidering more cornplex charge dependent modulation 
set-up (see e.g. Charge asymmetric modula- 

tion may also play a role explaining the discrepancy 
between the prediction of this model and low energy 
positron fraction PAMELA data (13j . 

A better fit can be obtained assuming a steeper 
spectrum at high energy, 70 — —2.5 rather than —2.42 
which worsen the fit of high energy Fermi-LAT data 
(see blue line in Fig. [1]) but may still be compati- 
ble with the current statistical and systematic errors. 
Forthcoming Fermi-LAT data and analysis will, there- 
fore, be crucial. Also this model, however, face serious 
problems when compared with low energy positron 
fraction PAMELA data. 

A consistent interpretation of all available exper- 
imental data in the framework of single component 



conventional CR propagation models is therefore still 
challenging. 



III. EXTRA COMPONENT 
INTERPRETATION 

Another viable interpretation of Fermi-LAT data 
consists in assuming the presence of two components 
in the CRE spectrum. This possibility was already 
risen in lOj and extensively discussed in [l3| and in a 
number of subsequent papers. In this scenario, an ex- 
tra electron and positron component of astrophysical 
(see e.g. [ll,[il]) or exotic origin (see e.g. [11 113,1111) 
is added to the conventional one. The spectrum of the 
extra-component must be significantly harder that the 
conventional one and have an exponential cutoff at 
Ecut — 1 TeV such to match high energy Fermi-LAT 
and H.E.S.S. data, while the low energy background 
electron should be softer than in the single component 
scenario. It was showed that if the extra-component 
is composed of electrons and positrons in compara- 
ble amount it allows to explain the rising behaviour 
of the positron fraction observed by PAMELA above 
10 GeV. In [l3] it was already noted as such scenario 
also allows to improve the fit of low energy pre- Fermi 
CRE data, especially those taken by AMS-01, respect 
to the single component conventional model. In that 
paper an injection spectral index —2.54 was adopted 
for the Galactic conventional component. 

The new Fermi-LAT low energy data favour a 
steeper spectrum. As shown in Figj^l the combination 
of Galactic electron background with source spectral 
index —2.7, modulated in force-field approximation 
with $ = 550 MV, and an electron -I- positron extra 
component with source spectral index P = — 1.5 expo- 
nentially cut at 1 TeV, provides an excellent fit of the 
whole CRE spectrum measured by Fermi-LAT as well 
low energy AMS-01 data and H.E.S.S. data above 1 
TeV. PAMELA positron fraction is also nicely repro- 
duced above 10 GeV while a small discrepancy re- 
mains below that energy (see insert in Figl2]) . In that 
figure the extra component source distribution has 
been assumed to coincide with that of Galactic pul- 
sars. We checked, however, that Fermi and PAMELA 
data can also be reproduced with a source distribu- 
tion which resemble that expected from dark-matter 
annihilation (though with a slightly different choice of 
the relevant parameters). 

A much better fit of low energy PAMELA data is 
obtained if a larger value of the diffusion coefficient 
spectral slope 6 is adopted. For example we found 
that if (5 ~ 0.5 (Kraichnan like diffusion) as suggested 
by a recent analysis [l^ of CREAM B/C measure- 
ments [20] and 70 ~ -2.6 the whole PAMELA and 
Fermi-LAT data set can consistently be fitted (the 
figures corresponding to that model will be showed 
elsewhere) . 
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FIG. 1: Fermi-LAT CRE data lioll. as well as several other experimental data sets, are compared to the e~ +e'^ spectrum 
modeled with GALPROP (see The gray band represents systematic errors on the CRE spectrum measured by 

Fermi-LAT. The spectra are obtained with injection indexes 70 = —2.42/2.5 (red/blue) above 4 GeV and 70 = —1.6 
below that energy. The diffusion coefficient spectral slope is 5 = 0.33 for both models. In the insert the positron fraction 
for the same models is compared with experimental data. Modulation in accounted in the force field approximation 
assuming a potential $ = 0.45 GV. 



As a specific realisation of such scenario we 
considered a model in which the extra electron 
and positron component is produced by galac- 
tic pulsars as illustrated in details in '14l|. Ob- 
served pulsars within 3 kpc from the Earth, as 
extracted from the ATNF radio pulsar catalogue 
fhttp : / / www. atnf . csiro .au/ research / pulsar/ps rcat / ), 
are accounted. We tested that adding more distant 
pulsars have negligible effects on our results. For 
each of those pulsars we use the spin-down luminosity 
given in the catalogue and randomly vary the rele- 
vant parameter in the following representative ranges: 
1 < Ecnt < 1-8 TeV, 10 < r?e± < 40 % (e± produc- 
tion efficiency), 5 < (At/10'' yr) < 10 (delay of 6=*= 
release respect to the pulsar birth) and 1.5 < F < 2.0. 
These ranges of parameter are compatible with our 
observational and theoretical knowledge of particle ac- 
celeration in PWNe. The electron-|-positron spectra 
expected for any of those realizations, as well as the 
predicted positron fraction, are shown in Fig[3]in the 
case 6 = 0.33. It evident that are several realizations 
allowing to match all electron data as well as high 
energy PAMELA positron fraction data. 



IV. CONCLUSIONS 

We reported on possible interpretations for the cos- 
mic ray electron-plus-positron (CRE) spectrum mea- 
sured by Fermi-LAT including preliminary new data 
extending down to £^ ~ 7 GeV 0, [l^ . We discussed 
the case of a single Galactic electron component, and 
a two-component scenario which adds to the conven- 
tional electron flux a primary electron and positron 
component with a harder spectrum, which may be ei- 
ther of astrophysical or exotic origin. 

In the single component scenario, Fermi-LAT elec- 
tron data can be approximatively be described if an 
injection spectral index between —2.4 and —2.5 is as- 
sumed (for S = 0.33). Those models, however, do not 
allow a perfect match of Fermi-LAT data either at low 
or at high energy. Furthermore they are in sharp ten- 
sion with the PAMELA data on the positron fraction. 

We therefore considered models adding an extra 
electron and positron component to the conventional 
electron Galactic sea. We showed that such models 
allow to improve considerably the fit of Fermi-LAT, 
H.E.S.S. and AMS-01 electon data as well as to ac- 
count for the rising behaviour of the positron fraction 
observed, most conclusively, by PAMELA. 
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FIG. 2: This model assumes the presence of an electron and positron (charge symmetric) extra component with spectral 
index V = —1.5 and exponential cutoff at i?cut — 1-2 TeV coming from a continuos distribution of sources in the Galactic 
disk (showed as a solid black line). The conventional electron component has an injection spectral index 70 = —2.7 
{5 — 0.33) (showed as a dotted black line). Both components are consistently propagated with GALPROP. Solid lines 
are modulated in the force field approximation with $ — 550 MV; dashed lines represent the LIS. 



Acknowledgments 

The Fermi LAT Collaboration acknowledges sup- 
port from a number of agencies and institutes for both 
development and the operation of the LAT as well as 
scientific data analysis. These include the NASA and 
the DOE United States, CEA/Irfu and IN2P3/CNRS 
in France, AST and INFN in Italy, MEXT, and the 
K. A. Wallenberg Foundation, the Swedish Research 



Council and the Swedish National Space Board in 
Sweden. Additional support from INAF in Italy 
for science analysis during the operations phase is 
also gratefully acknowledged. D.G. is supported by 
the Italian Space Agency under the contract AMS- 
02.ASI/AMS-02 n.I/035/07/0 and partially by Uni- 
verseNet EU Network under contract n. MRTN-CT- 
2006-035863. 



[1] T. Kobayashi, Y. Komori, K. Yoshida and 

J. Nishimura, ApJ 601 (2004) 340. 
[2] M. Aguilar et ai, Phys. Reports 366 (2002) 331. 
[3] A. W. Strong, I. V. Moskalenko and V. S. Ptuskin, 

Ann. Rev. Nucl. Part. Sci. 57 (2007) 285 
[4] F. A. Aharonian, A. M. Atoyan and H. J. Volk, A&A 

294 (1995) L41. 
[5] M. Pohl and J. A. Esposito, ApJ 507 (1998) 327. 
[6] J. Chang et al. [ATIC Collaboration], Nature 456 

(2008) 362. 

[7] F. Aharonian et al. [H.E.S.S. Collaboration], Phys. 

Rev. Lett. 101 (2008) 261104. 
[8] F. Aharonian et al. [H.E.S.S. Collaboration], 

arXiv:0905.0105 [astro-ph.HE]. 
[9] O. Adriani et al. [PAMELA Collaboration], Nature 



458 (2009b) 607. 
[10] A. A. Abdo et al. 2009 [Fermi collaboration], Phys. 

Rev. Lett. 102, (2009) 181101. 
[11] I. Moskalenko, A. Strong, Adv. Space Res. 27 (2001b) 

717. 

[12] M.S. Potgieter, Adv. Space Research 19, (1997) 883 
[13] H. Cast and S. Schael, Proc. of the 31th International 

Cosmic Rays Conference (ICRC 2009), Lodz. 
[14] D. Grasso et al. [Fermi Collaboration], Astropart. 

Phys. 32 (2009) 140 
[15] D. Hooper, P. Blasi and P. D. Serpico, JCAP 0901 

(2009) 025. 

[16] P. Blasi, Phys. Rev. Lett. 103 (2009) 051104 

arXiv:0903.2794 [astro-ph.HE]]. 
[17] L. Bergstrom, J. Edsjo and G. Zaharijas, Phys. Rev. 



eConf C091122 



2009 Fermi Symposium, Washington, D.C., Nov. 2-5 



5 



I A Kobayashi (1999) ■ ATIC-1,2 (2008) 

- ^ HEAT (2001 ) X PPB-BETS (2008) 

- o BETS (2001) ▼ HESS (2008-09) 




10' 10' 10' 

E (GeV) 
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is compared with experimental data. Each gray top line represents the sum of all pulsars for a particular combination 
of pulsar parameters. The dashed (pulsars only) and solid (pulsars + GCRE component) blue lines correspond to a 
representative choice among that set of possible realizations. The dot-dashed (purple) line represents the contribution 
of Monogem pulsar in that particular case. Note that merely for graphical reasons, here Fermi-LAT statistical and 
systematic errors are added in quadrature. In the insert the positron fraction for the same models is compared with 
experimental data. Solar modulation is accounted as done in FigH) 
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